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ABSTRACT 

We present an untriggered search for optical bursts with the ROTSE-I telephoto array. Obser- 
vations were taken which monitor an effective 256 square degree field continuously over 125 hours 
to mjiOTSE = 15.7. The uniquely large field, moderate limiting magnitude and fast cadence of 
'^lO minutes permits transient searches in a new region of sensitivity. Our search reveals no can- 
didate events. To quantify this result, we simulate potential optical bursts with peak magnitude, 

TOp, at t=10 s, which fade as / = ^t~^ j where at < 0. Simple estimates based on observational 

evidence indicate that a search of this sensitivity begins to probe the possible region occupied by 
GRB orphan afterglows. Our observing protocol and image sensitivity result in a broad region 
of high detection efficiency for light curves to the bright and slowly varying side of a boundary 
running from [a^jTOp] = [—2.0,6.0] to [—0.3,13.2]. Within this region, the integrated rate of 
brief optical bursts is less than 1.1 x 10"^ s~^ deg~^. At ~22 times the observed GRB rate 
from BATSE, this suggests a limit on < 5 where 6opt and 6^ are the optical and gamma-ray 
collimation angles, respectively. Several effects might explain the absence of optical bursts, and 
a search of the kind described here but more sensitive by about 4 magnitudes should offer a more 
definitive probe. 

Subject headings: gamma rays: bursts, observations - ISM: jets and outflows - stars: variables: other - 
methods: data analysis 



Motivation 

Gamma-ray bursts remain one of the great mys- 
teries in astrophysics. Although there have been 
some concrete measurements of the energetics of 
some bursts through redshift determination (eg. 
(Metzger et al. 1997), (Kulkarni et al. 1998)), there 
is little firm knowledge of how the energy is pro- 
duced. In fact, the total production is still un- 
certain by approximately 2 or 3 orders of magni- 
tude because of the unknown level of postulated 
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collimated jets. Additionally, despite recent ad- 
vances in multi-wavelength detection strategies, 
the total number of GRBs studied optically re- 
mains small. While various observations have 
placed the internal-external shock scenario on a 
relatively firm footing (eg. (Akerlof et al. 1999), 
(Sari and Piran 1999)) for some bursts, it is not 
verified for most. Bright optical bursts are the 
expected signature of reverse external shocks, 
but have been ruled out for several gamma-ray 
bursts (Akerlof et al. 2000), (Kehoe et al. 2001), 
(Park et al. 1999), which argues against a uniform 
behavior. 

Gamma-ray bursts are believed to emit syn- 
chrotron or inverse Compton radiation from ma- 
terial moving at ultra-relativistic velocities. The 
resultant strong Lorentz beaming of the emission 
will decrease as the shocked material slows down. 
If the ultra-relativistic bulk flow is a collimated 



jet, radiation at wavelengths longer than gamma 
rays, believed to be produced by this slower ma- 
terial, will be emitted through a larger solid an- 
gle (Rhoads 1998). This suggests a population of 
orphan optical bursts with timescales similar to 
GRBs, but more frequent and with no gamma-ray 
signature. 

With the unfortunate demise of the Compton 
Gamma-Ray Observatory, the rate of GRB de- 
tections overall has decreased dramatically. The 
physics of jets may prove valuable in the ongoing 
effort to observe multiwavelength emission from 
their progenitors and in a wider variety of condi- 
tions than studied so far. Conversely, the study 
of the relative rates of gamma-ray and "orphan" 
optical bursts may uniquely probe the open ques- 
tion of coUimation. Whether useful constraints 
can be obtained in this way is currently the sub- 
ject of some debate, some holding the view that 
strong limits are possible from such a measure- 
ment (eg. (Rhoads 1998)). On the other hand, 
(Dalai et al. 2002) have assorted these; efforts are 
unable to measure the coUimation angle. However, 
in (Dalai et al. 2002) only very late, slow cadence, 
and narrow-field searches were discussed. In addi- 
tion, the estimates of optical lightcurves are based 
heavily on an interpretation of a model which has 
yet to be confirmed for GRB's. 

Observational evidence for bursts at longer 
wavelengths than gamma-rays is at present very 
sketchy. One search for the expected higher 
rate of X-ray afterglows has been performed 
(Grindlay 1999) with negative results. On the 
other hand, mysterious optical transients have 
been observed in deep, narrow-field searches 
(eg. (Tyson et al. 2001)). Although selected 
based on spectral criteria rather than photomet- 
ric (lightcurve) critera as in our paper, the SNe 
optical search of (vandon Berk et al. 2001) has 
recently revealed a new type of AGN at 17th 
magnitude exhibiting characteristics reminiscent 
of GRB optical afterglows. This is close to the 
ROTSE-I detection threshold, and as a back- 
ground (Gal- Yam et al. 2002) will have to be un- 
derstood for any GRB orphans search. 

A substantial complication to these kinds 
of searches arises from the fact that the rates 
of rapidly varying optical transients are poorly 
understood. Study of eruptive variables with 
timescales of variation of order ~lh has been 



limited to searches with very non-uniform sky 
coverage. The catalog of known rapidly eclips- 
ing and pulsating systems has a similar limita- 
tion. The rates of burst-like AGN activity, such 
as from SDSS J124602.54-F011318.8, are also very 
poorly understood. Aside from being interesting 
in their own right, these transients are sources of 
background for GRB optical counterpart searches, 
whether triggered or untriggered. 

Method 

For all of these reasons, we have devised an op- 
tical burst search strategy based almost solely on 
considerations taken directly from data. Enough 
GRB late optical afterglows have been mea- 
sured to indicate that lightcurves of the form 

/ = , with at ~ —1, are typical. When 

extrapolated back to times nearer the burst phase, 
the observed afterglow lightcurve of GRB 970228, 
which is fairly typical of those observed so far, 
indicates emission brighter than 16th magnitude 
for the first 30 minutes. While the non-detections 
of optical bursts imply this is not the typical sce- 
nario (Kehoe et al. 2001), (Akerlof et al. 2000), 
(Park et al. 1999), the observed optical burst 
from GRB 990123 was actually brighter than this 
extrapolation would indicate. This sensitivity is 
attainable in 3 minute exposures with ROTSE-I 
(Kehoe et al. 2001), and it places the observation 
of an optical burst lightcurve brighter than this 
limit over 30 minutes just within the telescope's 
capability. We therefore take this lightcurve sig- 
nature as our fundamental search criterion. 

Because this search is limited by field area cov- 
ered and integrated observation time, it makes 
sense to maximize both in our study. The hy- 
pothesized 30 minute window of detectibility for 
ROTSE-I sets the number of different pointings we 
can accomplish while we expect an optical burst 
to be visible which, in turn, sets the cadence of 
our observations. We will discuss this more in 
the Observations section, but ROTSE-I can ac- 
complish 2 pointings in the allotted time. In gen- 
eral, for a rapidly fading lightcurve, the cadence 
of a wide-field search must increase with decreas- 
ing aperture. This permits fewer pointings by 
the smaller telescope but this is compensated for 
by the larger field of view. Despite the modest 
sensitivity of ROTSE-I, we are aided by its rela- 
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tively large etendue, D^il, a measure of sensitivity 
for such searches. This product varies relatively 
little over a wide range of instruments because 
the useful solid angle for fixed focal plane area 
decreases as as the aperture, D, increases. 
Thus, under the assumptions made in this paper, 
ROTSE-I is roughly equivalent to a 0.5m, f/1.8 
system. At the BATSE fluence sensitivity of 10"^ 
to 10"'^ ergs/ cm^, there are approximately two 
bursts per 4tt steradians per day. If the optical 
emission fills a cone 5 times wider than that for 
the gamma-rays, then a detectable optical burst 
potentially occurs once per 160 hours in a single 
ROTSE-I field of 256 square degrees. 

The uniquely large field and fast cadence of the 
observations taken for this analysis push transient 
searches in a very different part of observational 
phase space than have been pursued before. To 
both tune our search strategy and estimate our 
sensitivity to the sought after signature, we have 
employed a simple Monte Carlo simulation. This 
simulation has three main functions: (1) repro- 
duce the photometric behavior of ROTSE-I mea- 
surements, (2) explore our search efficiency for a 
range of input optical burst lightcurve parame- 
ters, and (3) quantify the search sensitivity given 
the observing protocols and limiting magnitudes 
of the actual data taken. We ignore in this sim- 
ulation model dependent assumptions about rela- 
tivistic beaming or spectral evolution so that our 
result can be interpreted as generally as possible. 
We will describe this simulation more in the Anal- 
ysis section. 

Observations 

Our search for optical bursts utilizes data taken 
with the ROTSE-I telephoto array (Kehoe at al. 2001), 
which consists of four telescopes arranged in a 
2x2 configuration, each with an 8 degree field of 
view and 14" pixels. The major challenge in our 
analysis is the successful rejection of ~1.5 x 10^ 
non-variable sources for over 2000 observations 
in each pointing, as well as suppression of astro- 
physical transients and instrumental backgrounds. 
At a bare minimum, any candidate for an optical 
burst must be well-detected in at least two obser- 
vations. This stipulation, in conjunction with the 
30 minute window constraint, largely dictates our 
observational approach. 



We have taken the data used in this analysis in 
two main pointing modes, stare and switch. The 
first mode consists of repeated images of a sin- 
gle field each night with no repointings to other 
locations. This stare data was acquired in two 
different periods. In mid-December of 1999, two 
separate fields were monitored on different days, 
while in mid-April, we covered another field ex- 
tensively, originally to monitor the X-ray nova 
XTE J1118+480 (Wren et al. 2001). Camera 'd' 
was not operational during half of these observing 
nights, and camera 'c' was inoperable on one. The 
total stare data set comprises an effective monitor- 
ing of 256 square degrees in 80 s exposures for 47.5 
hours. The observation dates and coordinates for 
these fields, and the total area covered, are given 
in Table 1. 

While the stare data provides excellent tem- 
poral coverage of any transients which may oc- 
cur, the field-of-view probed during our crucial 30 
minute window is unduly restricted. As mentioned 
above, the finite visibility window of the hypothe- 
sized optical bursts dictates that we maximize the 
area covered in the allotted time. With ROTSE-I, 
we can fit two pointings into this timeframe. We 
employ a protocol in which five 80 s exposures are 
taken at each pointing, followed by a repointing at 
a second location where another five exposures are 
taken. The telescope then returns to the first field 
and restarts. This leaves blocks of five 80 s expo- 
sures with 7 minute gaps for each of the two point- 
ings. One peculiarity of this switch mode occurs 
when we co-add pairs of images later - namely, 
each cycle has an un-coadded fifth image. This 
extra image alternates between preceding or fol- 
lowing the two CO- adds. We have chosen two of 
our standard sky patrol fields which pass through 
the zenith during the late summer and which have 
a galactic latitude greater then 20 degrees to avoid 
Galactic extinction and overcrowding. The data 
were taken from late August through early Oc- 
tober 2000 during periods chosen to avoid bright 
moonlight (see Table 1). 

For the observations presented in this analy- 
sis we instituted a dithering observing procedure 
to suppress backgrounds stationary on the CCD. 
This involved a Z-shaped rastering by about 2.5' 
on a side (~10 pixels) during the observing se- 
quence to relocate bad pixels in celestial coordi- 
nates in consecutive frames. 
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Data Reduction 

Reduction of the raw images taken for this 
study involves an initial calibration of the images, 
through a co-addition and recalibration phase, to 
lightcurve construction. This analysis generated 
14,000 raw images which were processed to ap- 
proximately 8000 epochs and 2.5 x 10^ photomet- 
ric measurements. The backgrounds and process- 
ing time incurred necessitated data reduction as 
the images were taken, while adding the ability to 
diagnose and, if possible, correct for various issues 
of image quality. 

The initial calibration of raw images is de- 
scribed in (Kehoe et al. 2001) and involves three 
main steps: correction of images, source extrac- 
tion, and astrometric and photometric calibration. 
Images are corrected through the subtraction of 
dark current and flat-fielding. The darks, flats 
and bad pixel lists are generated at the begin- 
ning of the three major periods of data collection. 
The corrected images are analyzed using SExtrac- 
tor (Bertin and Arnouts 1996) to produce object 
catalogs. To speed up processing, we have opti- 
mized the choice of clustering parameters while 
not significantly degrading the ability to extract 
sources in crowded environments. When an ob- 
ject aperture contains a bad pixel, we flag the ob- 
ject in that observation as bad. The photometric 
and astrometric calibration involves matching the 
resultant list of sources to the Hipparchos cata- 
log (H0g et al. 1998). As will be described in the 
Analysis section, our search depends on the sta- 
tistical study of lightcurve data. For our final cal- 
culated search efficiencies to be correct, we must 
reproduce the observed photometric fluctuations 
seen in actual ROTSE-I data. We find that an ir- 
reducible 0.01 mag fluctuation, in addition to the 
normal statistical one, reproduces the observed be- 
havior well. We assign this as the minimum sys- 
tematic error for all observations. 

Final Image Processing 

Once the astrometry for each frame is estab- 
lished, we can co-add images in pairs to improve 
our sensitivity. We do this instead of taking ex- 
posures of twice the length to increase the dy- 
namic range available in our 14-bit CCDs, to help 
suppress bad pixel and other instrumental back- 
grounds, and to permit a later augmentation of 



the search at twice the temporal granularity. For 
each field per night (ie. 4 for stare, 8 for switch), 
we independently number frames starting at '1' for 
the first image. We co-add each even numbered 
frame to its immediately preceding odd numbered 
frame if it is adjacent in time. The even numbered 
frame is mapped to the odd numbered frame using 
a bilinear interpolation of pixel intensities to best 
avoid position resolution degradation and image 
artifact creation. During the co-addition, intensi- 
ties of known bad pixels in the first and second 
images are replaced with the median of their 8 
neighbors and arc given zero weight. Pixels in 
the final co-added image with weight < 2.0 are 
rescaled to homogenize the image. Since the first 
and second images are dithered with respect to 
each other, bad pixel locations generally have cor- 
rect image data from one of the two images. As 
a further precaution, once the final co-added im- 
age is itself clustered and calibrated as described 
above, found objects are flagged if a bad pixel falls 
within their 5 pixel wide aperture. Typical image 
sensitivities of ttirotse = 15.7 were obtained for 
CO- added images. 

The next step of the processing chain involves 
removing images which exhibit instrumental prob- 
lems. Images must satisfy a Sct limiting magnitude 
of at least 14.75, regardless of whether they are co- 
adds. This generally removes epochs from early or 
late in a night which are affected by twilight, al- 
though it also rejects images taken when weather 
is not of good observation quality. The position 
resolution of an image must also be less than 0.15 
pixel (~2 arcsec). This removes images which did 
not CO- add well, due usually to some large obstruc- 
tion (eg. a tree late in the evening as the field 
moves towards the horizon) or rare slipping of the 
mount clutch which compromises the astrometric 
warp map. 

The last step before lightcurve construction in- 
volves the elimination from the data sample of 
epochs or whole nights exhibiting instrumental 
problems. For nights with TV operational cameras, 
we omit epochs when fewer than — 1 cameras 
are taking observation quality images. We require 
whole nights to have at least two contiguous hours 
of images that pass quality cuts, where contiguous 
means that no two consecutive epochs are missing. 
The surviving nights generally exhibit very stable 
quality over those observations we included. To 
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quantify our observations, wc calculate for each 
good night the median limiting magnitude of the 
images used, and the total observing time. These 
are indicated in Figure 1. 

Lightcurve Construction 

For each camera, the good object lists for each 

night arc matched to provide lightcurvcs for all 
detected sources. We remove any object which is 
not seen in at least one consecutive pair of obser- 
vations, while we retain all isolated observations 
of objects passing this selection. To avoid suspect 
photometry, intermittent sources arising from de- 
blending fluctuations, or moving objects, we flag a 
source's observations when the measured position 
is more than one half pixel from the mean source 
location. 

After we have created a filtered list of calibrated 
lightcurves, we must perform an internal correc- 
tion to remove systematic photometric mismea- 
surements due to such difficulties as the presence 
of thin haze over a portion of a frame, or increased 
vignetting due to a sticking shutter. More serious 
problems occur when power-lines or trees cross an 
image or when shutters obstruct a portion of an 
image. In each of these cases, we determine either 
a small relative photometric correction and sys- 
tematic error, or in more extreme situations flags 
to indicate the problem and exclude the image 
area. We begin by selecting a set of good template 
sources in each matched list. For each source to 
enter this list, we count those observations which 
have good photometry. If these criteria are satis- 
fied for more than 75% of the epochs in which the 
source location was imaged, the object becomes a 
template source and we calculate its median mag- 
nitude in good observations. Each image is then 
divided into 100 pixel square subfiles, and pho- 
tometric offsets for each template source having 
statistical error less than 0.1 mag in these subtiles 
are calculated. The number of template sources, 
rii, within a typical subtile, i, is around 50. A 
relative photometry map is then calculated to be 
the median of these offsets, o^, for each subtile. 
The standard deviation of the offsets, c7o,i, which 
for good subtiles is typically around 0.03 magni- 
tudes, is used to calculate a systematic error 
on Oj. In addition, subtiles are classified as bad if 
rii < 5, or if (jo,i > 0.1 magnitude. These selec- 
tions remove obstructions and large photometric 



gradients due to out-of-focus foreground objects, 
which constituted approximately 2% of the data 
sample. Each observation of each object is then 
corrected by the Oj using a bilinear interpolation, 
and the Si is added in quadrature to that obser- 
vation's systematic error. If the subtile is flagged 
as bad, then the object's observations are likewise 
flagged. 

The lightcurve analysis proceeds from those ob- 
servations of a particular source which pass basic 
quality cuts. The source must neither be satu- 
rated, nor near an image boundary. The aperture 
magnitude must be properly measured. The po- 
sition of the source in good observations must be 
less than 0.5 pixels from the source's mean coor- 
dinates. The inefficiency incurred for these selec- 
tions is negligible. There must also be no known 
bad pixels within the source aperture, and the 
subtile containing the source must have a well- 
measured photometric offset. The former ineffi- 
ciency is approximately 10%, while the latter is 
2% as stated earlier. 

Analysis 

Lightcurve Selection 

Based on our observationally motivated model, 
we look for lightcurves which may have timescales 
for the brightest optical emission which is similar 
to the gamma-ray emission of a GRB. In addi- 
tion, we may expect a decaying lightcurve with 
at < 0. This signature is unlike any other known 
type of optical transient, and our lightcurve se- 
lection exploits this. We begin by cutting on sim- 
ple statistical variables designed to efficiently iden- 
tify a fast-rising transient with a power-law decay, 
while removing backgrounds. These backgrounds 
come in two kinds. The first originate from astro- 
physical processes such as low amplitude pulsat- 
ing and eclipsing systems, irregular fast variables, 
flare stars, and asteroids. The second kind is domi- 
nated by various instrumental transients like those 
resulting from unregistered bad pixels, poor pho- 
tometry, or moderate mount clutch failures not 
identified by our image quality cuts. 

In order to obtain an a priori selection, we 
first studied the instrumental background by per- 
forming a grid search on the April 9th stare data 
of camera 'b' (2% of the total search sample) in 
terms of the maximum variation amplitude. A, 
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and the significance of that variation, cta- We 
took as comparison the 13 echpsing and pulsating 
systems we observe in this sample down to varia- 
tion amplitudes of approximately 0.1 magnitude. 
We explored the region of A = [0.1. .0.25] and 
CTA = [2.0. .5.0], and we selected cuts of A > 0.1 
and (Ta > 5.0 which minimized the total back- 
ground while retaining substantial efficiency for 
real, low amplitude variables. 

These variables generally come in two main 
types: pulsating stars such as the five indi- 
cated in Figure 2, and eclipsing systems such 
as the three illustrated in Figure 3. Of these 
8 variables, only one is previously identified as 
such: the Algol binary BS UMa (ie. ROTSEl 
J112541. 62+423448.8) (Skiff 1999). The rest are 
low amplitude variables difficult to identify in pho- 
tographic plates. ROTSEl Jl 12037.63+392100.4 
and ROTSEl J113536.74+384557.4 may be a 5 
Scu and an RRc, respectively, but we have not 
fully classified the variables in this sample at this 
time. However, averaging over four fields, we esti- 
mate the rate for variables with amplitude > 0.1 
mag and period < Id to be of order O.ldeg"^ 
for eclipsing systems and 0.2deg~^ for pulsating 
systems. Several variables with longer timescale 
variation (ie. 0.7 — 3d) are also observed, and 
the approximate rate of these is 0.2deg~^. Note 
that the data is out of the galactic plane. 

We further investigated the instrumental back- 
ground behavior in switch data taken in July, 
2000. These fields are not included in the search 
due to a substantial worsening of the problem with 
the tracking of the mount. This was later al- 
most completely fixed, greatly reducing the back- 
ground in subsequent switch data. The remain- 
ing backgrounds largely consisted of very stable 
sources with one anomalous photometric measure- 
ment. Cutting on a calculated from obser- 
vations not including the most significant varia- 
tion, Xci '^^^ effective against these. A candidate 
flare-star, ROTSEl J160542.74+350016.3, was ob- 
served in outburst twice in this July data as shown 
in Figure 4. 

We examined the effect of this selection on our 
hypothesized optical burst signal. We simulated 
the expected GRB signature with a simple Monte 
Carlo generating a variety of lightcurves with no 
emission during the first 10 s, and then peaking at 
nip from 6.0 to 15.0. A power-law fading with at 



from -0.05 through -3.0 followed thereafter. A cut 
of Xc > 3.0 is very efficient for the simulated bursts 
and removes most backgrounds passing our other 
cuts. Note that the calculation of Xc requires an 
object be seen in at least three good observations. 
This in turn means our search will be insensitive 
to variations lasting less than about 7 minutes. 
All of the simulated optical bursts passing our Xc 
cut also had a A > 0.5 magnitude. We tighten 
our selection accordingly to remove most of the re- 
maining low amplitude periodic variables seen in 
the data background samples. Thus our lightcurve 
cuts require an overall significant variation with 
A > 0.5 and cta > 5.0, and substantial variabil- 
ity in the rest of the lightcurve with Xc > 3.0. 
Although designed with a simple power-law time 
decay in mind, these cuts are efficient for a wide 
range of optical burst lightcurves. 

At this stage, our candidate sample was domi- 
nated by three major sources of background. The 
first consisted largely of bright, moderate am- 
plitude pulsating and eclipsing systems, most of 
which have not yet been cataloged. The second 
category consisted of non- varying sources present 
in most observations but poorly measured in more 
than one observation. For the most part, these 
latter are not flare stars. These categories of can- 
didate are easily measured on other nights and 
so rejected. Given the results of our simulated 
lightcurves, these are highly unlikely to correspond 
to GRB optical counterparts. The last category 
arises from the fact that the outer region of each 
camera field has relative photometry corrections 
which start to break down where interpolation be- 
comes extrapolation. We discard these regions in 
our search, and we thereby incur a 2.5% ineffi- 
ciency. 

Figure 5 indicates, for a co-add sensitivity of 
15.7 magnitude, the 50% efficiency contour for our 
simulated bursts after our final lightcurve cuts. 
GRB 970228 and GRB 990123 are shown for com- 
parison. We see that the selection is efficient for 
simulated bursts in a particular subset of the nip 
vs. at parameter space. The boundary of this 
region is relatively sharp, and moves vertically as 
a night's typical image sensitivities increase. The 
increased sensitivity to steep lightcurves with the 
stare protocol results from the uninterrupted na- 
ture of the observing. 
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Search Results 

After the selection mentioned in the previous 
section, 50 candidate Ughtcurves remain. Twenty 
six of these candidates are located along a CCD 
column in camera 'c' with improper charge trans- 
fer and so rejected. Hand-scanning and compari- 
son to images without candidates revealed that 12 
of the remaining 24 lightcurves were due to bad 
pixels. On three nights (000824, 000903, 000906), 
there was an intermittent 15% increase in the 
number of hot pixels observed in data. The reason 
for this behavior is not understood. Although it 
affects a small portion of image real-estate, these 
pixels were not properly dark-subtracted, flagged 
or removed in co-addition. One set of candidates 
happens when two bad pixels are near enough that 
they land on the same a, 6 in consecutive frames. 
The other candidates consist of a hot pixel land- 
ing between two nearby bright stars, overlapping 
a star or stars just below reliable detection, or fol- 
lowing a satellite trail in the next image. 

Ten of the remaining 12 candidates are within 
15" (~1 pixel) of objects in the USNO A2.0 cat- 
alog (USNO 1998) brighter than 17th mag. in R. 
By rejecting these, we incur an inefficiency dom- 
inated by the loss of the image region contained 
in ~ 6 pixel area around all of the USNO source 
positions. Since the largest occupancy per camera 
is about 50,000 stars, this cut incurs a 7% ineffi- 
ciency. 

The last two candidates clearly move slowly 
in our images. These match the expected coor- 
dinates of the asteroids (1719) Jens and (17274) 
2000 LC16, and we reject them. Thus, we find no 
optical bursts in this study. 

Our data sample was taken at several different 
sensitivities and both stare and switch pointing 
modes, and this will alter the sensitivity of our 
search from that shown in Figure 5. Using the 
simple Monte Carlo already described, we gener- 
ated over 20 million bursts folding in the fractions 
with co-add sensitivities and pointing mode given 
in Figure 1 . A region of [at , nip] parameter space 
is efficiently accepted which lies on the bright and 
long side of a rough line between [-2.0,6.0] and 
[-0.3,13.2]. This result is shown in Figure 6. A fi- 
nal additional inefficiency of 20% for our selection 
comes from three main sources: bad pixel removal 
(e = 0.9), bad image regions (e = 0.98), removal of 



previously known sources (e = 0.93), and removal 
of field edges (e = 0.975). 

Discussion 

We have performed the most extensive untrig- 
gered search for optical bursts in wide-field data to 
date. After study of 1331.4 deg^ days, we find no 
candidates down to a typical limiting sensitivity of 
rriROTSE = 15.7. This sensitivity is of the order 
necessary to study the potential rates for optical 
bursts from GRBs based on simple beaming hy- 
potheses. At a maximum efficiency of 80%, we ac- 
cept transients which are brighter and longer than 
a boundary running from [Q;t,mp] = [—2.0,6.0] 
to [—0.3,13.2]. In this region, we therefore rule 
out an integrated optical burst rate greater than 
1.1 X 10~^ s~^ deg~^. At a rate ten times greater, 
we reject bursts bounded by [at,mp] = [—2.6,6.0] 
to [-0.3,14.3]. 

As there is no empirical information about the 
optical properties of orphan afterglows, any in- 
terpretation of this result in the context of GRB 
beaming is uncertain. If we consider that the op- 
tical emission is coUimatcd into some jet angle, 
9opt, which may be different than the gamma- 
ray collimation angle, 6j, then the relative rate of 
gamma-ray to optical bursts may indicate some- 
thing about the relative magnitudes of these an- 
gles. Compared to the observed GRB rate from 
BATSE of 5 X 10^1° s^^ dcg"^ our limit on opti- 
cal burst rates translates into a limit on < 5_ 

This result seems to conflict with the much 
higher rates predicted in the analysis of (Frail et al. 2001). 
There may be several reasons for this. We have 
assumed that all of the bursts observed by BATSE 
would have produced optical bursts of sufficient 
intensity to be observed above 16th magnitude. 
The relationship of optical emission during this 
early phase to gamma-ray emission is not known 
well enough to rigorously justify this assumption. 
In this case, however, a search of the kind de- 
scribed in this analysis but deeper by 2-4 magni- 
tudes should provide a more incisive probe, given 
the distribution of GRB fluences relative to those 
we have discussed in this paper. 

Although most GRBs with both X-ray and ra- 
dio afterglows seem to have optical counterparts, 
it is still possible that obscuration at the source 
is dimming a large number of the optical bursts 
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wc would otherwise see. For instance, there is ev- 
idence for as much as 5 mag extinction for GRB 
970828 (Groot et al. 1998). Ultimately, the ques- 
tion of extinction at the source is still tied up with 
the identity of the progenitor. In coUapsar mod- 
els (Woosley 1993), substantial extinctions could 
result. 

Lastly, we have ignored the specifics of how 
the changing Lorentz factor of the burst ejecta 
affects the very early optical emission. This be- 
havior is complex, and very model dependent, but 
for very off-axis viewing the observed optical flux 
is likely to peak later, and at a lower value (eg. 
(Granot et al. 2002)), than might be suggested by 
extrapolating back from the late afterglow. Ob- 
servationally, many optical afterglows have exhib- 
ited breaks in their lightcurves which may indi- 
cate analogous behavior for nearly on-axis viewing 
(Frail et al. 2001). In a couple of cases, the early 
lightcurve exhibits a time decay index of ~ —0.8 
((Stanek et al. 1999), (Stanek et al. 2001)) which 
puts them in a relatively sensitive region of Fig- 
ure 6. Any optical burst with such a slope and 
brighter than about 12th mag at it's peak would 
have been observed, and bursts with peak bright- 
nesses of ~ 13th mag would be detected with lower 
efficiency. As one moves more off-axis, the slope 
at early times becomes shallower until eventually 
it may rise to the break point. The suppression of 
this early emission will result in a decreased abil- 
ity to observe these bursts. However, it should be 
kept in mind that, regardless of when the peak 
occurs relative to the initial event, our search will 
be sensitive to these lightcurves if the emission is 
brighter than 16th mag in a way similar to the 
optical bursts we have simulated (ie. detection in 
at least three images with overal A > 0.5). In 
any case, a limiting magnitude deeper by ~ 4 mag 
would substantially improve the chances for detec- 
tion, especially in light of observed GRB optical 
afterglow brightnesses. 

This search technique can be improved and ex- 
tended in several ways. The main source of back- 
ground was instrumental. To consider much larger 
numbers of observations, it is important to reduce 
the number of bad pixels and shutter malfunc- 
tions, and avoid visual obstructions. Another im- 
provement would be to increase the sophistication 
of the relative photometry calibration. Most im- 
portantly, the next step for an untriggered optical 



burst search would be to consider deeper fields, 
such as available with the ROTSE-III 0.5m tele- 
scopes. This kind of search will likely be looking 
for more gradual optical variation than sought in 
this paper. For robust identification, it may be 
necessary to obtain spectral verification of candi- 
date orphans which are more likely to be found at 
the sensitivity of ROTSE-III. The search must be 
made immune to the kind of background engen- 
dered in SDSS J124602.54-F01 1318.8, in particu- 
lar. 

This analysis has implications for our ability to 
extend our GRB studies. We have already used 
the techniques described here to search for optical 
counterparts for several archival BATSE triggered 
bursts (Kehoe et al. 2001). In addition, they are 
a crucial step towards our goal of near real-time 
detection with the ROTSE-III telescopes the opti- 
cal transients associated with GRB triggers from 
satellite-based gamma-ray observations. 

ROTSE is supported by NASA under SR&T 

grant NAG5-5101, the NSF imdcr grants AST 97- 
03282 and AST 99-70818, the Research Corpora- 
tion, the University of Michigan, and the Plane- 
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ported by the Department of Energy under con- 
tract W-7405-ENG-36. Work performed at LLNL 
is supported by the Department of Energy under 
contract W-7405-ENG-48. 
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Fig. 1. — Distribution of integrated observing 
times vs. limiting magnitude for the 64 square 
degree field of a single camera. The limiting mag- 
nitudes range from 15.0 to 16.1 for the stare data, 
and 15.4 to 16.1 for the switch data. The total 
observing time for this size field is 20.8 days. 
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Fig. 2. — Five single-night lightcurves for candi- 
date pulsating variables from the April 2000 stare 
data. None of these are previously identified as 
variable. Observations indicated are good obser- 
vations after relative photometry. Errors are stat. 
+ sys. 
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Fig. 3. — Three lightcurves for candidate eclipsing 
binary variables from the April 2000 stare data. 
Only ROTSEl J112541. 62+423448. 8 was previ- 
ously known to be variable. Observations indi- 
cated are good observations after relative photom- 
etry. Errors are stat. + sys. 
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Fig. 4. — A candidate flare star from the July 
2000 test switch data. Observations indicated are 
good observations after relative photometry. Er- 
rors are stat. + sys. Two outbursts were observed 
4 days apart with amplitudes of approximately 0.5 
mag and durations of about 30 min. The rest of 
the lightcurve, including the intermediate day July 
5, show the source at a constant 12.18 m. The 
lightcurves are arbitrarily aligned to peak at the 
same time. 



13 




16 I I 

-1 -2 -3 

Power-law Index 

Fig. 5. — Sensitive region for observing proto- 
cols. The red outline delineates 50% efficiency for 
switch data, while the green outline indicates the 
stare data 50% efficiency. The efficiency turn-on 
for a particular limiting magnitude is quite steep. 
Both are calculated for co-add image sensitivies 
of 15.7. The increased efficiency in stare data for 
steeply decaying lightcurves is indicated. The lo- 
cations of GRB 970228, extrapolated from after- 
glow data (Galama et al. 2000), and GRB 990123 
(Akerlof et al. 1999) are shown for comparison. 
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Fig. 6. — Cumulative sensitivity to simulated 
sources of given nip and at- Contours are gen- 
erated by folding the observing time as given in 
Figure 1 with the protocol efficiency contours such 
as given in Figure 5. Contours giving the 95%, 
75%, 30% and 10% efficiency levels are specified. 
The gradient in efficiency is primarily due to the 
varying protocol and limiting magnitudes of the 
ROTSE-I observations. The efficiencies are given 
relative to a rate of 1.1 x 10~* s~^ deg~^. 
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year and month 


dates 


a 




f2(dec;2) 


dwell-time x area 


pointing mode 


Dec, 1999 


16 


4h 


+ 15° 


768 


4.355 X lO^deg^s 


stare 




17 


5h24m 


-hll.4° 


512 




stare 


Apr, 2000 


9,10,13,14 


llh42m 


+44° 


768 


3.9375 X lO^'deg^s 


stare 




15-17 






(1024) 






Aug, 2000 


24,25,31 


Oh 


+30° 


1024 


4.2799 X lO^deg^s 


switch 


Sep, 2000 


1-3,5,6,28,29 












Oct, 2000 


1-4,6 












Aug, 2000 


24,25,31 


23h 


+15° 


1024 


2.8488 X lO^deg^s 


switch 


Sep, 2000 


1-3,5,6,29 












Oct, 2000 


1-4 













Table 1: Itemization of observation properties. 

The coordinates of fields and total area covered 
are given, as well as the dates and observing pro- 
tocol for each. 
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